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The London penetration depth is a fundamental quantity of a superconductor. It plays a key role in determining 
and uncovering many properties of a superconductor, including homogeneity, anisotropy, isotope effects, importance 
of quantum and thermal fluctuations, and facets of the nature of superconductivity in a particular material. In 
recent years experimental data on the temperature, dopant concentration, magnetic field and oxygen isotope mass 
dependence of the penetration depth became available for a variety of cuprate superconductors. Here we analyze and 
discuss the experimental data, guided by the generic phase diagram of the cuprates, depicted in Fig.l. After passing 
the so called underdoped limit {pu ~ 0.05), where p is the hole concentration, Tc reaches its maximum value T™ at 
Pm ~ 0.16. With further increase of p, Tc decreases and finally vanishes in the overdoped limit po ~ 0.27 [1,2]. There is 
the line Tc (jp) of finite temperature phase transitions, separating the superconducting and non-superconducting states, 
with critical endpoints at pu and Po ■ Here Tc vanishes and the cuprates undergo at zero temperature doping {p) tuned 
quantum phase transitions. As their nature is concerned, resistivity measurements reveal a quantum superconductor 
to insulator (QSI) transition in the underdoped limit and in the overdoped limit a quantum superconductor 
to normal state (QSN) transition. Another essential experimental fact is the doping dependence of the anisotropy. 
In tetragonal cuprates it is defined as the ratio 7 = Xc/^ab of the London penetration depths due to supercurrents 
flowing perpendicular (Ac ) and parallel {Xab ) to the ab-planes. Approaching the QSN transition 7 remains finite, 
while at the QSI transition it tends to infinity [3,4]. When 7 remains finite the system exhibits anisotropic but 
genuine three dimensional (3D)-, while 7 ^ 00 implies 2D-behavior. The resulting competition between anisotropy 
and superconductivity raises serious doubts whether 2D mechanisms and models, corresponding to the limit 7 = 00, 
can explain the essential observations of superconductivity in the cuprates [4] . There is mounting evidence that close 
to the phase transition line Tc (p) thermal fluctuations dominate, while quantum fluctuations dominate both, the 
QSI and QSN transitions. Furthermore, due to the 3D to 2D crossover, tuned by the rise of 7 with reduced dopant 
concentration, these fluctuations are enhanced. For these reasons, mean-field treatments, including the BCS theory 
are expected to apply far from the critical line only. This singles out the low temperature region around optimum 
doping p^. 
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FIG. 1. Schematic phase diagram of cuprate superconductors. Variation of and 7 (T = 0) = Ac (T = 0) /\ab (T = 0) with 
hole concentration p. Quantum phase transitions occur at the endpoints p = pu — 0.05 and p = po — 0.27 of the critical line 
Tc (p). At pu a two dimensional quantum superconductor to insulator (2D-QSI)- and at po ^ 0.27 3D quantum superconductor 
to normal state (3D-QSN)- transition occurs. 

In the mean-field approximation the London penetration depth of an anisotropic superconductor in the Meissner 
state is given by [5,6] 



A| (T) TT^Jki^ J \wf\ \ J A, dEk /B^ 



dSpyVPi and |vf| arc respectively the surface element of the Fermi surface, Fermi velocity in direction i and the 
magnitude of the Fermi velocity. / (Ek) is the Fermi function and Ak the energy gap in direction k. The i index refers 
to the principal crystallographic directions a, b and c. The second term which is negative, describes the decrease of 
1/A| (T) caused by the thermal population of Bogoliubov quasi-particlc levels with energy Ek, and it is this quantity 
where the anisotropy and magnitude of the energy gap enters. In this approximation 1/Xf (T) vanishes close to Tc 
as 1/Af (T) = 1/A|) (1 — T/TcY with v = 1/2, while their is mounting evidence for v ~ 2/3 in the experimentally 
accessible regime. Thus close to Tc where thermal fluctuations dominate the mean field treatment fails. Otherwise, e.g. 
at sufficiently low temperatures and far away from the QSI and QSN transitions the neglect of fluctuations appears 
to be justified. Here Eq.(l) reduces for noninteracting quasiparticle excitations around the four d-wave nodes, which 
dominate the leading low temperature behavior, to [7] 



The Fermi velocities vp and V2 enter the quasiparticles excitation energies Ei~ = /i-y/v^ifc^"+v^fcf and refer to velocities 
along directions normal and tangential to the Fermi surface at each node, d is the mean interlayer spacing along 
the c-axis. The velocity ratio vf/v2 is a fundamental material parameter which measures the anisotropy of the 
quasiparticle excitation spectrum. This scenario is not restricted to the penetration depth. It predicts simple power- 
law temperature dependencies in the thermodynamic and transport properties at sufficiently low temperatures. For 
example, the penetration depth measurements find that 1/A^^ exhibits in the clean limit and at low temperature a 
linear temperature dependence [8,9], in agreement with Eq.(2). The NAIR relaxation rate exhibits the expected 
temperature dependence [10]. The predicted effect of impurities in giving rise to a universal thermal conductivity 
[11,12] has been confirmed [13]. The clean-limit specific heat varying as appears to have been observed [14 16]. 

For a spherical Fermi surface, one recovers for the zero temperature penetration depth the standard result, 
1/A^ (T = 0) = 47rne^/ (mc^), where n is the number density of the electrons in the normal state [17]. In recent 
years the penetration depth has been the subject of intense experimental investigation in high- Tc compounds [18 24]. 
It has become the practice to associate the anisotropy with anisotropic effective masses m* and to interpret the 
experimental data in terms of the London formula 

A|(0) m*c2 ^-"^ 
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by introducing the number density of the superfiuid. In a real superconductor, the ionic potential modifies the 
spherical Fermi surface of free electrons drastically. In this case it is not evident what should be taken for the efi'ective 
mass m* and the number density rig. Although, using Eqs.(2) and (3), we can define the ratio ns/m*, which often 
has been used to interpret experimental results, in terms of 

1 f^^^v^^Wn^^ 



A?(0) n^hc'^ J ^ |VF 

This relation shows that ns/nii is just a way of parameterizing experimental results, with no discernible connection to 
the band mass or carrier concentration. Indeed in the mean-field approximation 1 /A? (0) is determined by normal state 
properties, namely the integral of the Fermi velocity over the Fermi surface. Noting that in high- Tc superconductors 
band structure calculations [25] and ARPES [26,27] studies uncovered drastic deviations from the free electron picture 
it is evident that interpretations based on the ratio ns/m* obscure the origin of the anisotropy 7jj = and the 

doping dependence of the zero temperature penetration depth, as well as the isotope effect on this quantity. Indeed, an 
inspection of Eq.(2) leads to the conclusion that the anisotropy stems from flat portions of the Fermi surface, while the 
doping dependence reflects that of the Fermi surface. On the other hand , there is little doubt about the importance of 
residual electron-electron and electron-phonon interactions, not accounted for in Eq.(2). Quantifying these interactions 
is difficult in the normal state of the cuprates, given the lack of well-defined single-particle excitations as revealed by 
various experiments. Contrariwise, well-defined quasiparticle excitations do exist in the superconducting state, and a 
description of the low temperature state in terms of superfluid Fermi liquid theory is believed to apply. Fermi-liquid 
corrections account for the Fcrmi-liquid interactions between electrons. In the superconductor their effect is the 
renormalization of the Fermi velocity ratio in terms of Vj?/v2 o^YpjNi [7], where is the bare value entering 

Eq.(2). The comparison oid\-'^/dT at T = in Bi2212, evaluated with Eq.(2) and the ARPES estimate for vf/v2 [26] 
with the value deduced from penetration depth measurements [28] points to a substantial Fermi liquid renormalization, 
namely « 0.3, due to interactions between the nodal quasiparticles in the superconducting state. To explore the 
doping dependence of the renormalization we invoke the empirical relation Tc d/dT (A^^ (0) /A^;, (T)) |^ ^ « —0.6 [9], 
which applies to variety of cuprates with Tc ranging from 30 to 130K and dopant concentrations extending from the 
underdoped to the optimally doped regime. It implies with the empirical relation, Tc oc 1// A^;, (0) [23] , established 
for underdoped cuprates, and Eq.(2) that o^nfI^i is nearly doping independent. 

Much less attention has been devoted to the renormalization of the Fermi velocity due to electron-phonon inter- 
action. Isotope substitution, i.e. the exchange of ^^O by is a suitable probe, whereby the lattice parameters 
remain essentially unaffected [29,30], while the phonon frequencies associated with the mass of the oxygen ions, or 
more generally, the oxygen lattice degrees of freedom are modified [31]. In recent years the isotope effect on the 
zero temperature penetration depth has been investigated, using a variety of techniques. In Table I we listed the 
experimental estimates for the relative change 

upon isotope exchange for various cuprates and MgB2, where n =^^0, m =^®0 in the cuprates and n =^^B, m =^°B 

in MgB2. The data also reveals that within experimental accuracy AA^^ [t^ /AJ;^ {t^ = AA^^ (0) /A^^ (0). At zero 

temperature and taking the anisotropy of the quoted materials into account (Ac >> Aa « Ab) the mean-field expression 

(2) reduces to 

' dS^F^,. (6) 



(T = 0) 7^2nc2^/2 

Since the lattice parameters remain essentially unaffected [29,30] by isotope exchange, while the dynamics, associated 
with the mass of the respective ions, are modified, the substantial isotope effect on the zero temperature penetration 
depth requires a renormalization of the normal state Fermi velocity v^? Vi? = Vi?/(l + /) where v^? is the bare 
velocity and / the electron-phonon coupling constant which changes upon oxygen isotope exchange in the cuprates 
but remains nearly unaffected by boron isotope substitution in MgB2. However, in the Migdal-Eliashberg [33] (ME) 
treatment of the electron-phonon interaction the coupling constant / is independent of the ionic masses and assumed 
to be small [34,35]. This is true if the parameter ujof /Ef is small, where ujq is the relevant phonon frequency and 
Ef the Fermi energy. Thus the isotope effect on the penetration depth is expected to be small, of the order of the 
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adiabatic parameter 7 = ojq/ Ep << 1. The ME theory retains terms only of order 0. Cuprates, however, have Fermi 
energies much smaller than those of conventional metals [36] so that 7 is no longer negligible small. Thus the large 
oxygen isotope effects on the zero temperature in-plane penetration depth in the cuprates, listed in Table I, poses a 
fundamental challenge to this understanding and calls for a theory that goes beyond ME [37,38] . 





Tc{K) 


T 


AA^, (t) (t) 


Ref 


YBaaCuaOr-a 


89.0 


10 


0.05 


[20] 


Y0.7 Pro.3Ba2Cu307_5 


60.6 


5 


0.13 




Y0.6 Pro.4Ba2Cu307_5 


45.3 


5 


0.11 




YBaaCuaOy.i (film) 


89.3 


4 


0.05 


[21] 


La2-£!;Sra;Cu04_(-5, x = 0.08 


19.5 





0.10(2) 


[19] 


La2-a;Sra;Cu04+5, a; = 0.086 


22.4 





0.08(1) 




Bii.6Pbo.4Sr2Ca2Cu30io+5 


107 





0.05 


[18] 


MgB2 


38.5 





0.02(2) 


[24] 



m 



Table I: Experimental estimates for AA^^ (t) (t) = ("A^^ (t) -™ A^^ (t)) /"A^^ (t) with n =^^0, 
=^^0 in the cuprates and n =^^B, m =^°B in MgB2 



Indeed, the relative shifts AA^^ {t^ /A^j, {t^ are substantial and surprisingly close to (^^Mo-^^Mq) /^^M^ = 
0.125. This differs fundamentally from the behavior of optic O - phonon frequencies. The expected behavior, uj oc 

M~^/^, was confirmed in YBa2Cu307_5 by measuring the frequency shift of the transverse optic phonons (copper- 

1/2 

oxygen stretching modes), yielding Aw/w « —0.06, in agreement with Aw/w « (^^Mq/^^Mq) — 1 = —0.057 [31]. In 
any case the observed oxygen isotope effect on the zero temperature penetration depth uncovers together with Eri.(6), a 
substantial renormalization of the normal state Fermi velocity due to oxygen lattice degrees of freedom, while in MgB2 
this renormalization due to boron isotope exchange turns out to be marginal. This renormalization is also expected 
to affect the superconducting properties. Taking the empirical relation Tc d/dT (A^^ (0) /A^;, (T)) ~ —0.6 [9] for 
granted, Eq.(2) implies that 



AT. 



AA^ao) 



(7) 



where vp/v2 is the bare ratio, renormalized with respect to electron- phonon coupling. Noting that close to optimum 
doping ATc/Tc is negligible small, the oxygen isotope exchange uncovers a substantial electron-phonon renormal- 
ization of vp/v2, characterizing the quasiparticles in the superconducting state. Indeed, close to optimum doping 

A (vf/v2^ I (vfM) = - AA^^ (0) /A^b (0) holds and this quantity varies from 0.05 to 0.11 (see Table I). This effect, 

providing direct evidence of the relevance of electron-phonon coupling in the superconducting state, should be ob- 
servable with ARPES. Noting that thermally excited quasiparticles destroy superconductivity by driving l/A^^ (T) to 
zero, we can estimate by extrapolating Eq.(2) to 1/A^j, (T) = 0. This yields T^A ss 1 and confirms Eq.(7). In this 
context it is interesting to note that details of the Fermi-surface topology of deuterated k-(BEDTTTF)2Cu(NCS)2 
have been measured as a function of pressure and compared with equivalent measurements of the undeuterated salt. 
The data suggest that the negative isotope effect observed on deuteration is due to small differences in Fermi-surface 
topology caused by the isotopic substitution [32]. 

Close to the phase transition line Tc (p) thermal critical fluctuations, neglected in mean-field treatments, dominate 
the thermodynamic properties. Approaching the phase transition line around p = from below (see Fig.l), there 
is mounting evidence that the critical behavior of homogeneous cuprates falls in the experimentally accessible tem- 
perature regime into the 3D-XY universality class [4,39-52]. Here critical 3D-XY fluctuations dominate because the 
fluctuations of the vector potential are strongly suppressed due to the small value of the effective charge of the pairs 
[39]. In the 3D-XY universality class the transition temperatTirc^ T,;, the critical amplitude of the specific heat A~ and 
the critical amplitudes of the penetration depths A^o arc universally related by [50] 



167r3 



\2 \2 \2 A- 



(8) 
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R~ ~ 0.815 is a universal number. The critical amplitudes are defined as c = {A~/a)t~°' , t = 1 — T/Tc and 
A? = )\ot~'^, where a and i' are the critical exponents. Although Tg, A~ and A?o depend on the dopant concentration, 
isotope exchange etc., universality implies that this combination does not. Hence this relation puts a crucial constraint 
on the microscopic theory of superconductivity in cupratcs. To illustrate this point it is instructive to consider the 
doping dependence of Tg, A^;, ~ AaoAbo 7o = ^co/^abo and A~ close to the 2D-QSI transition (see Fig.l). Here Tg, 
1 / A^fjQ and 70 are known to scale as oc 1 /Xlbo ^ /To oc p — p„ so that the critical amplitude of the specific heat 
singularity vanishes according to Eq.(8) as A^ cx (p — Pu) oc T^. This is consistent with the specific heat data for 
underdoped La2-a:Sra;Cu04 and underdoped Tl2Ba2Cu06+5 [53]. Furthermore, the universal relation implies that 
the relative changes upon isotope substitution satisfy the relation [54] 

3^1- 3 ^ A? ■ ^ ' 

i=a,b,c »0 

It explains why close to the transition temperature and optimum doping (p = pm) a substantial isotope effect on the 
penetration depths is compatible with a negligible effect on Tc- 

Approaching the 2D-XY-QSI transition in the underdoped regime (see Fig.l) a crossover to the universal relation 

[50] 

T = ^ (10) 

takes place and the cuprates correspond to a stack of independent sheets of thickness ds- R2 is a universal number. The 
flow to 2D-XY-QSI behavior is experimentally well confirmed in terms of Uemura's plot [23] . It is a characteristic 2D 
property and also applies to the onset of superfluidity in "^Hc films adsorbed on disordered substrates where it is well 
confirmed [55]. Although T^, ds and A^f, (0) depend on dopant concentration, isotope substitution, etc., universality 
implies that this relation does not. This puts yet another constraint on the microscopic theory of superconductivity 
in the cuprates. Furthermore it yields for the relative changes upon isotope exchange the relation 



AT, _ Ms AA^,, (0) 
To ds A2,(0) 



(11) 



applicable close to the imdcrdopcd limit [p = p^)- Although the experimental data arc rather sparse in the imdcrdoped 
regime [54,56], suggestive evidence for an isotope effect on the effective thickness dg of the superconducting sheet 
emerges, namely Ads/ds = {^^ds —^^ dg) /^^dg ~ 0.03 [54]. Estimates for dg can be derived from the crossing point 
phenomenon in the temperature dependence of the magnetization for various applied magnetic fields, applied parallel 
to the c-axis. In a 3D anisotropic superconductor, falling into the 3D-XY universality class, the magnetization data 

1 /2 

plotted in terms rUz/Hz versus T will cross at T,, where z is along the c-axis. Consistency with this behavior was 

found in a variety of cupratcs [50]. On the contrary, in a 2D superconductor, corresponding to a slab of thickness 
dg, the crossing point occurs in the plot rriz versus T at the Kosterlitz-Thouless transition temperature Tkt, where 
ruz oc —ksTKr/ {^odg). Even though bulk cuprates are strictly 2D only close to the 2D-QSI transition at p = p„ 
(sec Fig.l), in the highly anisotropic materials, such as Bi-2212 and Tl-1223, 2D crossing point feature have been 
observed and used to estimate ds , yielding values close to the c-axis lattice constant [50] . Since the lattice parameters 
remain essentially unaffected [29,30] by isotope exchange, while ds does, the substantial relative change Adg/ds = 
(^^ds —^^ dg) /^^dg ~ 0.03 [54] imcovcrs again the relevance of electron-lattice coupling. Going further, by combining 
Eqs.(7) and (11), extrapolated to the underdoped regime, we obtain the approximate relation 



Ads ^ [^Fh2 



ds 



(12) 



which provides additional evidence for the coupling between superconducting properties and lattice degrees of freedom. 
The approximate nature of this relation stems from the fact that the d-wave quasiparticle scenario does not hold down 
to the 2D-QSI transition due to the neglect of quantum fluctuations, associated with the phase of the order parameter. 
Indeed, the linear-in-T temperature dependence of 1/A^^ simply follows from the existence of a 2D-QSI transition. 
The result is, {Xab (0) /\ab {T)) — 1 oc T/Tc, holds for d- and s-wave pairing [9]. 

Additional evidence for electron-lattice coupling emerges from the combined oxygen isotope and finite size effects. 
Due to inhomogeneities of extrinsic or intrinsic origin, cuprates are homogeneous over a finite extent only. Thus a 



5 



finite size effect [57] is expected to occur, whereby the correlation volume cannot grow beyond the volume of the 
homogeneous domains. When 3D-XY critical fluctuations dominate there is the universal relationship [50] 



Ai (T)A, (T) 



(13) 



between the London penetration depths A; and transverse correlation lengths 5* in directions i and j. In the presence 
of inhomogeneities with length scales the ^| = 1^1"", where t = T/Tc — 1, cannot diverge but are bounded by 



(14) 



A characteristic feature of the resulting finite size effect is the occinrcnce of an inflection point at Tp,^ in 1/Aj (T) \j (T) 
below Tc, the transition temperature of the homogeneous system. Here 



and Eq.(13) reduces to 



Ww^kBTp^ 1 



A, (T)A, (T) 



Lk 



(15) 



(16) 



In the homogeneous case 1/ (Aj (T) Aj (T)) decreases continuously with increasing temperature and vanishes at 
while in the presence of inhomogeneities it exhibits an inflection point at Tp^ < Tc, so that 



d 



Xi {T)Xj{T) 



/dT 







(17) 



Since the experimental data for the temperature dependence of the penetration depths is available in the form Xah 
and Ac only, we rewrite Eq.(16) as 



Ww^kBTp^ {Xab {T) Ac {T))t=t^^^ 

Lnh = ^ 



(18) 



Apart from the inflection point, an essential characteristic of a finite size effect is the finite size scaling function [58]. 
In the present case and for Xab it is defined in terms of 



Xpab 
Xab (T) 



\t\ =9c{y), y = sign {t) \t\ 



^Oab 



sign {t) 



t 



(19) 



For t = T/Tc — 1 small and Lc oo, so that ±y — > oo, it should tend to gc {y — oo) = 1 and gc {y oo) = 0, 
respectively, while for t = and Lc 7^ it diverges as 



9c (y ^ 0) = goc \y\ " = goc 



t 



(20) 



whereby (XQab/Xab (Tc, L))'^ = gQc\tp,.\'^ = .^oc^oabZ-^c- As expected, a sharp superconductor to normal state 
transition requires domains of infinite extent. Moreover at tp^, y = 1 and there is an inflection point because 
d{XQab/Xab{T,L)f' /dt = 0. Since the scaling function gc{y) depends on the type of confining geometry and on 
the conditions imposed (or not, in the case of free boundaries) at the boundaries of the domains, this applies to 
the amplitude g^c as well. In Fig. 2a we displayed the microwave surface impedance data for A^^ (T = 0) /A^^ (T) 
and d (A^^ {T = 0) /A^j, (T)) /dT versus T of a high-quahty Bi2Sr2CaCu208+A- single crystal taken from Jacobs et 
al. [59]. The solid curve indicates the leading 3D-XY critical behavior of the homogeneous system, while the data 
uncovers a rounded transition which occurs smoothly. This behavior, together with the occurrence of an inflection 
point around Tp^ « 87K, where c^(A^j {T ~ 0) /A^j, (T)) /dT exhibits an extremum, points to a finite size effect. 
With Xab {T = 0) = I8OOA obtained from fiSR measurements [60], Tp^ « 87K and Xl^ (T = 0) /A^^ (TpJ = 0.066 
we obtain with the aid of Eq.(18) the estimate Lc « 68A. Although the spatial extent of the homogeneous domains 
along the c-axis appears to be of nanoscale only, the small critical amplitude of the transverse correlation length. 
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^t^^ = Lc (1 - TpjTcf'^ « 2.2A, makes the 3D-XY critical regime (A^^ (0) /\\ (T) cx (1 - T/T,f'^) attainable. An 
additional and essential characteristic of a finite size effect appearing in the temperature dependence of the in-plane 
penetration depth is the consistency of {Xoab/^ab (T))"^ \t\~'^ versus t/ \tpj with the shape and limiting behavior of 
the finite size scaling fimction (sec Eq.(19)). In Fig. 2b wc displayed (Xoab/Xab (T)) \t\~'^ versus t/ \tpj. The apparent 
agreement with the aforementioned characteristic behavior of this function, provides strong evidence for a finite size 
effect, due to the limited extent Lc of homogeneous superconducting domains along the c-axis. Clearly, such a finite 
size scaling analysis, performed on one set of data for one particular sample and material only, cannot distinguish 
between an intrinsic or extrinsic origin of the inhomogeneity. Noting that this behavior was found in a variety of 
cuprates and for data obtained with different techniques [61], one is lead to the conclusion that inhomogeneities, giving 
rise to a finite size effect, are yet another facet of the mystery of superconductivity in the cuprates. Clearly this finite 
size effect is not restricted to the penetration depth but should be visible in other thermodynamic properties. In the 
specific heat it leads to a rounding of the peak and its consistency with a finite size effect was established for the data 
taken on YBa2Cu307_5 high quality single crystals [50]. In these samples the domain size was found to range from 
300 to 400 A. Furthermore nanoscale spatial variations in the electronic characteristics have also been observed in 
underdoped Bi2Sr2CaCu208+5 with scanning tunnelling microscopy (STM) [62-65]. They reveal a spatial segregation 
of the electronic structure into 3nm diameter superconducting domains in an electronically distinct background. 




T ( K ) y = " tp 

FIG. 2. (a)Microwave surface impedance data for A^j (T = 0) /Xlt, (T) (Q) and d (A^j (T = 0) /A^f, (T)) /dT (□) versus T of a 
high-quality Bi2Sr2CaCu2 08+i single crystal taken from Jacobs et al. [59]. The solid line is A^f, (0) /Xl^ (T) = 1.2 (1 - T jTcf'^ 
and the dash-dot line its derivative with Tc = 87. 5K, indicating the leading critical behavior of the homogeneous system. The 
dotted line is the tangent to the inflection point at Tp « 87K, where d (A^t (0) /Aq^ (T)) jdT is maximum; (b) Finite size scaling 
function g (y) = {Xoab/Xab {T))^ \t\~" versus y = t/ \tp\ for the data shown in Fig.2a. The solid line in the inset is Eq.(19) with 
goc = 1-6. 



Supposing that the limiting length scales change upon isotope exchange, we obtain from Eq.(18) the relation 

which matches Eq.(ll), applicable in the 2D limit where the limiting length is set by ds, the thickness of the in- 
dependent sheets. Furthermore, this relation opens the possibility to probe the coupling between superfluidity and 
lattice degrees of freedom close to criticality, where mean-field treatments fail. Indeed, given the fact that the lat- 
tice parameters remain essentially unaffected [29,30] by isotope exchange, a purely electronic mechanism requires 
ALc = 0. The effect of oxygen isotope siibstitiition on the inhomogeneity induced finite size effect has been explored 
in Yi_3;Pra;Ba2Cu307_5 [66]. From the resulting estimates, listed in Table II, several observations emerge. First, 
Lc increases systematically with reduced Tp^. Second, Lc grows with increasing x and upon isotope exchange (^^O, 
-'^^O). Third, the relative shift of Tp^ is very small. This reflects the fact that the change of Lc is essentially due to 
the superfiuid, probed in terms of 1/A^^. Accordingly, ALc/Lc « AA^j {TpJ /A^^ {TpJ for a; = 0, 0.2 and 0.3. 
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X 





0.2 


0.3 




-0.000(2) 


-0.015(3) 


-0.021(5) 




0.12(5) 


0.13(6) 


0.16(5) 




0.11(5) 


0.15(6) 


0.15(5) 




89.0(1) 


67.0(1) 


52.1(2) 




89.0(1) 


66.0(2) 


51.0(2) 




9.7(4) 


14.2(7) 


19.5(8) 


^«ic (A) 


10.9(4) 


16.0(7) 


22.6(9) 



Table II: Finite size estimates iovATpjTp^, AL^/L^ AKb (^Pc) 1Kb (^Pc)> ^'^^'pc ^^^'pc ^^L^ and ^^L^ for an i^O 
content of 89% taken from [66]. 

To appreciate the implications of these estimates, we note again that for fixed Pr concentration the lattice parameters 
remain essentially unaffected [29,30]. Accordingly, an electronic mechanism, without coupling to lattice degrees of 
freedom implies ALc = 0. On the contrary, a significant change of L,. upon oxygen exchange uncovers the coupling 
between the superfluid, probed by A^;,, and the oxygen lattice degrees of freedom. A glance to Table II shows that the 
relative change of the superconducting domains along the c-axis upon oxygen isotope exchange is significant, ranging 
from 12 to 16%, while the relative change of the inflection point T.p^ is an order of magnitude smaller. For this reason 
the significant relative change of at fixed Pr concentration is accompanied by essentially the same relative change 
of A^f,, which probes the superfluid. This uncovers unambiguously the existence and relevance of the coupling between 
the superfluid and oxygen lattice degrees of freedom. Furthermore, this behavior agrees with the isotope effect on 
ds [54], the limiting length scale close to the 2D-QSI transition. Potential candidates for the relevant lattice degrees 
of freedom are the Cu-0 bond-stretching-type phonons showing temperature dependence, which parallels that of the 
superconductive order parameter [67]. 

An additional probe to unravel the mystery of superconductivity in the cuprates is the response to a magnetic 
field. In the early discussion of the symmetry of the order parameter, Yip and Sauls [68] proposed that the angular 
position of the gap nodes could be probed by a measurement of the magnetic field dependence of the penetration 
depth. In the local limit and for T^O, they predicted the hnear relationship, (A (if = 0, T = 0) /A (if, 0))^ — 1 oc —if, 
where the factor of proportionality is independent of temperature. Several experimental groups tried to verify this 
prediction, but failed to identify a linear if term which scaled with temperature according to the theory [69-72]. 
On the other hand, calculations based on a d-wave model, treated in the quasi classical approximation, suggest that 
(A (if = 0, T = 0) /A (if, 0))^ — 1 oc —^/H [74]. However, in these treatments fluctuations have been neglected. Close 
to the phase transition line Tdp), where thermal fluctuations dominate, the in- plane penetration depth scales as 
Kb ^ ic = iabh^ while the magnetic field applied parallel to the c-axis scales as ifg oc ^'o/Cob- Thus the in-plane 
penetration satisfies at Tc the scaling form 

fW„ (22) 



A^JT^ff,)^^"- 

revealing that a superconductor is dramatically influenced by an applied magnetic fleld. This behavior can be un- 
derstood by noting that in an applied magnetic field the correlation length cannot grow indefinitely. For nonzero 
magnetic field ifg there is the limiting length scale Lh^ — \/^o/ {o-Hc) with a ~ 3.12 [61], related to the average 
distance between vortex lines. Indeed, as the magnetic field increases the density of vortex lines becomes greater, but 
this cannot continue indefinitely. The limit is roughly set on the proximity of vortex lines by the overlapping of their 
cores. Due to this limiting length scale the phase transition of a homogeneous superconductor is rounded and occurs 
smoothly. At T = and close to the 2D-QSI transition the in-plane penetration depth scales as A^;, oc where z is 
the dynamic critical exponent of this transition. With ifc oc ^o/^ab yields the scaling form 



A (0,0) 



1 oc -H^/"^. (23) 



^A(ffc,T = 0), 

Contrariwise at T = and close to the 3D-QSN transition the in-plane penetration depth scales as A^^ a ^^^^ so that 



A (0,0) 

A(ffe,T = 0) 



- 1 oc -ff^^+1)/^ (24) 
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Since the order parameter is assumed to be a complex scalar, these scaling forms hold for both, s-wave and d-wave 
pairing. Taking the evidence for a 2D-QSI transition with z = 1 and a 3D-QSN transition with z = 2 into account 
[4], the scaling form {X {H = 0,T = 0) / X {H, 0))^ — 1 oc —sfH'c is expected to hold in the underdoped regime, while in 
the overdoped Hmit (A (i7 = 0, T = 0) /A 0))^ - 1 cx -i/c should apply. In any case more extended experimental 
investigations are required, including samples covering the full doping range, to overcome the present impasse. Another 
property suite to shed light on the critical properties of the quantum transitions is the magnetic field dependence of 
the zero temperature spcicific heat coefficient. At T = and close to the 2D-QSI or 3D-QSN transitions it scales as 
(c/T)j,^Q oc oc Hc^^^^^^. The data taken on La2_a;Sra;Cu04 [75] points to {D — .z) /2 ~ 1/2, irrespective of the 

dopant concentration. This suggests respectively, 2; = 1 for the 2D-QSI and z = 2 for the 3D-QSN transition. 

The doping dependence of the zero temperature penetration depths provides another link between the quantum criti- 
cal behavior in the underdoped and overdoped regimes. In Fig.3 we displayed 1/A^j, (0) versus p ~ x for La2-xSra;Cu04 
taken from Panagopoulos et al. [76]. Close to the quantum phase transitions 1/A^j, (0) scales as l/A^^ (0) cx 
[50], where 5 — p — pu at the 2D-QSI and S = po — p &t the 3D-QSN transition. The solid line indicates the crossover 
from a 2D-QSI transition with z = 17 = 1 to a 3D-QSN transition with z = 2 and V = 1/2. While the flow to 
the 2D-QSI transition is apparent, the data does not extend sufficiently close to 3D-QSN criticality to confirm this 
crossover unambiguously. In any case, it emerges that the properties of the ground state are controlled by the crossover 
from the 2D-QSI to the 3D-QSN critical point. For this reason it can be understood that away from these quantum 
critical points, around optimum doping {p = Pm, sec Fig.l), quantum fluctuations are suppressed to the extent that 
Bogoliubov quasi-particle features are observable and mean-field treatments represent a reasonable starting point. 
A 2D-QSI transition with z = 1 and F = 1 coincides with the theoretical prediction for a 2D disordered bosonic 
system with long-range Coulomb interactions [77-79] . A potential candidate for 3D-QSN criticality is the disordered 
d-wave superconductor to disordered metal transition at weak coupling considered by Herbut [80] [38], with z = 2 
and V = 1/2. energy. Here the disorder destroys superconductivity, while at the 2D-QSI transition it localizes the 
pairs. 



3.0 




FIG. 3. l/Xat, (0) versus p ~ x for La2-a;Sra;Cu04. •: experimental data taken from Panagopoulos et al. [76] The solid line 
indicate the crossover from a 2D-QSI transition with z =T' = 1 io & 3D-QSN transition with z = 2 and V = 1/2. 

To siimmarizc, in the regime where mean-field treatments are expected to apply, the substantial isotope effect on 
the zero temperature penetration depth, established by a variety of experimental techniques, implies a renormalization 
of the normal state Fermi velocity due to a electron-lattice coupling, beyond the band velocity and the ME theory. 
This coupling also affects the superconducting properties and should be observable with ARPES. Close to the critical 
line Tc {p) and the critical endpoints, either thermal or quantum fluctuations, not included in mean-fleld treatments, 
dominate. In these regions of the phase diagram the theory of quantum and thermal critical phenomena applies. Given 
the universality class of the respective transition, there are universal relations between critical properties, putting 
stringent constraints on the microscopic theory of superconductivity in the cuprates. Along the phase transition line 
there is in the experimentally accessible temperature regime mounting evidence for 3D-XY universality and for a 
3D-2D- crossover as the underdoped limit is approached. Here a 2D-QSI transition occurs. This crossover, measured 
in terms of the ratio 7 = Xc/Xab (see Fig.l) is well documented [3,4] and implies that superconductivity in the 
cuprates is a germinc 3D phenomenon. The resulting competition between anisotropy and superconductivity raises 
serious doubts whether 2D mechanisms and models, corresponding to the limit 7 = 00, can explain the essential 
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observations of superconductivity in the cuprates. Indeed, as the dopant concentration is increased, the cuprates 
undergo in the ground state a crossover from 2D-QSI to 3D-QSN criticahty. For this reason the observation of 
Bogoliubov quasi-particle features far away from these quantum critical points, around optimum doping {p = pm, see 
Fig.l), can be understood. Here quantum fluctuations do not dominate. Yet another facet of the mystery emerges 
from the evidence for a finite size effect in the temperature dependence of the in-plane penetration depth. Although 
the limiting length scales may depend on the history of the sample, their dependence on oxygen isotope substitution 
reveals and confirms the coupling between superconducting properties and lattice degrees of freedom. Although the 
majority opinion on the mechanism of superconductivity in the cuprates is that it occurs via a purely electronic 
mechanism, and lattice degrees of freedom are supposed to be irrelevant, we have shown that the oxygen isotope 
effect on the in-plane penetration depth uncovers yet another facet, the hitherto ignored coupling between lattice 
degrees of freedom and both, normal state and superconducting state properties, the existence of inhomogeneities 
giving rise to a finite size effect. Finally we note that a variety of other properties also display pronounced phonon and 
electron-lattice effects: superconductivity-induced lattice changes [45,51,73,81,82], superconductivity-induced phonon 
renormalization [67,83-87], tunnelling phonon structures [88,89], etc., give additional evidence of significant electron- 
lattice coupling. Furthermore we have seen that the occurrence of power law terms in the low temperature and low 
magnetic field dependence of the in-plane penetration depth do not necessarily single out d-wave pairing, but stem 
from fluctuations at work, associated with a complex scalar order parameter which is compatible with both, s- and 
d-wave pairing. 
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